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Abstract 

FerB from Paracoccus denitrificans is a soluble cytoplasmic flavoprotein that accepts redox equivalents from NADH or 
NADPH and transfers them to various acceptors such as quinones, ferric complexes and chromate. The crystal structure and 
small-angle X-ray scattering measurements in solution reported here reveal a head-to-tail dimer with two flavin 
mononucleotide groups bound at the opposite sides of the subunit interface. The dimers tend to self-associate to a 
tetrameric form at higher protein concentrations. Amino acid residues important for the binding of FMN and NADH and for 
the catalytic activity are identified and verified by site-directed mutagenesis. In particular, we show that Glu77 anchors a 
conserved water molecule in close proximity to the 02 of FMN, with the probable role of facilitating flavin reduction. 
Hydride transfer is shown to occur from the A-pro-S position of NADH to the solvent-accessible si side of the flavin ring. 
When using deuterated NADH, this process exhibits a kinetic isotope effect of about 6 just as does the NADH-dependent 
quinone reductase activity of FerB; the first, reductive half-reaction of flavin cofactor is thus rate-limiting. Replacing the 
bulky Arg95 in the vicinity of the active site with alanine substantially enhances the activity towards external flavins that 
obeys the standard bi-bi ping-pong reaction mechanism. The new evidence for a cryptic flavin reductase activity of FerB 
justifies the previous inclusion of this enzyme in the protein family of NADPH-dependent FMN reductases. 
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Introduction 

There is an expanding list of bacterial flavoenzymes tliat display 
redox activities against diverse xenobiotics, such as heavy metal 
ions like chromate [I], nitro compounds [2], azo dyes [3] and 
carbonyl compounds including quinones [4]. These enzymes 
occur as homodimers or homotetramers with a non-covalently 
bound flavin mononucleotide (FMN) cofactor and display a bi-bi 
ping-pong mechanism with NAD(P)H as the presumed physiolog- 
ical reductant. Despite a progress in elucidating high-resolution 
structures of some representatives and understanding the under- 
lying flavin redox chemistry, key questions remain still open, in 
particular as regards the identity of natural substrates and the true 
significance in bacterial metabolism. 

The flavoprotein FerB (ferric reductase B, a product of gene 
4071) was originally isolated from the cytoplasm of the soil 
bacterium Paracoccus denitrificans as one of two major enzymes able 
to reduce Fe(III)-ligand complexes when NADH was the electron 
donor [5]. The protein is also active as a chromate reductase [5,6] 
and, to a substantially greater extent, as a quinone reductase [5,7]. 
Among the soluble quinone substrates tested, 2,3-dimethoxy-5- 
methyl- 1 ,4-bezoquinone (ubiquinone-0, UQ;0) yielded the highest 
tat/ K,„ specificity ratio, suggesting that FerB may, in vivo, interact 



with benzoquinone ring of UQ;10, a lipophilic electron carrier 
abundantly present in the cytoplasmic membrane of the bacteri- 
um. 

An initial Pfam domain analysis [7] placed FerB within the 
protein family of NADPH-dependent FMN reductases (Pfam 
database accession number PF03358). Some members of this 
family contain by themselves a non-covalently bound flavin 
cofactor (as does FerB) and are in fact quinone reductases, 
catalyzing two-electron reduction of quinones to respective 
hydroquinones [4] . Since modulation of the intracellular activity 
of quinone reductases in many instances affects cell responsiveness 
to redox stressors, these enzymes are believed to serve as protective 
devices for preventing semiquinone formation and subsequent free 
radical production or for quinone-mediated quenching of H2O2 
[8]. However, additional roles beyond detoxification are not 
excluded as well. For example, there is evidence that eukaryotic 
quinone reductases associate with the 20S proteasome in redox 
state dependent manner and control in this way proteasomal 
degradation of some transcription factors [9] . From the practical 
point of view the substrate specificity towards metal ions makes 
these enzymes good candidates for biodegraders of highly toxic 
inorganic contaminants such as chromate or uranyl [10]. They 
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may also find therapeutic application in prodrug processing by 
reductive enzymatic activation [11]. 

In this study, our main goal is to identily the amino acid residues 
of FerB responsible for flavin cofactor binding and/ or catalysis of 
its redox conversions. For this, the crystal structure of holoenzyme 
has been solved and the effect of mutating the strategically selected 
residues has been analyzed by means of fluorimetric titration of 
free flavin by apoenzyme, enzyme kinetic assays and spectro- 
potentiometry. Although there are no direct structural data on the 
enzyme-NADH complex, novel pieces of information about the 
interactions of FerB with NADH are provided by the observed 
stereochemistry of hydride abstraction from NADH and computer 
docking analyses. We also show that a simple steric modulation of 
the active site enhances its abiUty to reduce external flavin, which 
is otherwise negligible for the wild- type enzyme. 

Materials and Methods 

Gene Cloning, Protein Expression and Purification 

The cloning, expression, and purification of the FerB protein 
and its selenomethionine derivative with a C -terminal six-histidine 
tag (FerBHise) in E. colt strain BL21(DE3) pLysS cells (Invitrogen) 
were carried out as described in the paper [12]. 

Site-directed Mutagenesis 

Site-directed mutagenesis was carried out as specified in the 
QuikChange II Site-Directed Mutagenesis Kit (Stratagene). The 
pET21-ferB plasmid sensed as a PGR template; the primers are 
hsted in Table 81. After the mutagc'n(;sis protocol, the sequences of 
the new ferB constructs were verified by sequence analysis. The 
mutated plasmids were purified using the Plasmid Mini DNA 
Purification Kit (Qiagen) and transformed into the host expression 
strain E. coli BL21(DE3) pLysS (Invitrogen). Expression of the 
recombinant protein mutants and the purification procedure by 
Ni-IDA (GE Healthcare) chromatography were performed as 
described for the wHd-recombinant type enzyme [13]. Correct 
folding of all mutant enzymes was confirmed by circular dichroism 
spectroscopy. 

Crystallization, X-ray Data Collection and Structure 
Determination 

Crystallization, X-ray diflraction data collection and reduction 
of the native FerBHisR and its selenomethionine derivative are 
described in [12]. The crystal structure of FerBHisg protein was 
determined by multiple wavelength anomalous dispersion using 
selenium anomalous signal. Excellent data quality and resolution 
allow us to try to solve phase problem of FerBHisg using the 
advanced 3W-MAD protocol of Auto-Rickshaw: the EMBL- 
Hamburg automated crystal structure determination platform 
[14]. The Auto-Rickshaw structure determinative resulted in 
almost complete model of Se-Met FerBHise homodimer contain- 
ing 349 residues (out of 380) divided into 9 chains with Rl Rj„j, = 
0.2144/0.2682. Detailed description of the Auto-Rickshaw strat- 
egy in this particular case is described in [15]. Missing fragments 
were manually built using Coot [16] and this model of FerBHisg 
was re-refined against 1.40 A diffraction data and rebuild using 
PDB-REDO algorithm [17]. Final model of FerB (PDB ID: 
3U7R) contains 364 residues with R/Rj,^ = 0.16065/0.17386. 

Structure Alignments 

A search for sequence and structural analogy was done using 
Basic Logical Alignment Search Tool (BLAST) and DALI server 
[18]. Multiple alignment analysis and phylogenetic tree construc- 
tion was done using the PHYLIP program. 



SAXS (Small-angle X-ray scattering) 

The SAXS data were coUected on the BioS.AKS-1000, Rigaku 
at CEITEG (Brno, Czech Republic). Data were collected at 
290.15 K with X-ray beam wavelength X= 1.54 A. Sample to 
detector (PILATUS lOOK, Dectris Ltd.) distance was 0.4 m 
covering a scattering vector range from 0.008 to 0.65 A For 
buffer and sample one two-dimensional image was coUected with 
an exposure time of 30 min per image. The sample was measured 
at three concentrations: 5.0, 2.5 and 1.25 mg/mL. The buffer was 
identical to the buffer of the last step of the protein purification. 
Data reduction and the buffer substraction were performed using 
SAXSLab, Rigaku. Subtracted data were normalized to the 
protein concentration and truncated to 0.2 A~'. Evaluation of the 
solution scattering of the FerBHiss crystal structure and the fitting 
to experimental data was performed by CRYSOL [19]. The 
oligomeric state of the protein was evaluated using OLIGOMER 
[20] where the form-factors files were created by FFMAKER 
using dimeric and tetrameric assembly of the FerB crystal 
structure, both programs from ATSAS package [2 1] SAXS. Data 
collection and scattering-derived parameters are summarized in 
Table S2. 

Preparation of Apoflavoproteins 

Approximately 2 mg of purified FerBHis(, or its mutants in the 
buffer containing 50 mM sodium phosphate, 300 mM NaCl and 
10 mM imidazole, pH 8.0 were applied onto a 5 mL pre- 
equUibrated Ni Sepharose 6 Fast Flow column (GE Healthcare). 
After loading the protein, the column was washed 20 mL of 
starting buffer. Protein-bound flavins were then removed by 
washing the column with 100 mL of the same buffer additionally 
with 4 M KBr and 4 M urea (pH 8.0). The apo-form of the 
protein (apoFerBHisg) was eluted by the buffer containing 50 mM 
sodium phosphate, 300 mM NaCl and 500 mM imidazole, 
pH 8.0 (bufier C). 

Analytical Size Exclusion Chromatography 

0.1 mL of sample (0.5 mg/mL of proteins) was loaded onto 
Superose 12 (GE Healthcare) analytical and preparative column 
(30 cmxl cm) equilibrated by 25 mM Tris-HCl (pH 7.4). The 
flow rate after sample application was 0.5 mL/min. Marker 

proteins were albumin (67 kDa), ovalbumin (43 kDa), chymotryp- 
sinogen A (25 kDa) and ribonuclease A (13.7 kDa). 

DSC (Differential scanning calorimetry) 

Calorimetric measurements were performed on a VP-DSC 
MicroCalorimeter (MicroCal) with the cell volume 0.5 mL. Data 
were collected from 25 to 95 °C at a heating rate of 60 °C per 
hour. Protein concentrations were adjustixl to 2.5 mg/mL for 
FerBHisg and 4.3 mg/mL for apoFerBHisg in buffer C. AU 
solutions were degassed before measurements. Reference baseline 
was obtained by buffer vs. buffer scan and subtracted from the 
measured data. The obtained data were analyzed using the Origin 
software (Microcal Software). 

Dissociation Constants 

Dissociation constants of apoenzyme-flavin complexes were 
determined from fluorescence titration with excitation at 450 nm 
by adding of constant amount of the wild-type and mutant variants 
of apoFerBHisg into 25 mM Tris-HCl bufier (pH 7.4) with 
100 nM FMN or riboflavin. Decreasing emission intensities at 
519 nm were recorded on Luminescence Spectrotometer LS-50B 
(Perkin-Ekner). Series of different concentrations of FMN standard 
was used for setting of a calibration curve on fluorescence 
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intensity. The dissociation constants Kj were calculated according 
to [22] and converted to the binding Gibbs energy changes AGj - 
-RTln {\/Kjj (where R, gas constant, T, absolute temperature), 
from which the variation of Gibbs energy due to mutation (AAGy 
was computed as the difference in binding energy between mutant 
and wild-type enzyme (AAGj = AGj(mut) - AG4(wt)). 

Steady-State Kinetics 

The initial velocity of NADH oxidation was measured 
spectrophotometricaUy at 340 nm (s = 6.22 mM"' cm"') after 
the addition of variable concentrations of NADH to the solution 
containing the enzyme and a fixed concentration of electron 
acceptor (UQ-O or flavin) in 25 mM Tris-HCl (pH 7.4). The 
kinetic parameters ± standard error were determined by fitting 
the data to the appropriate kinetic equation using nonlinear 
regression and the program Microcal Origin (version 5.0). 

Determination of the Standard Redox Potential 

The standard redox potentials of FerBHisg and its mutants were 
determined at 25 °C according to the method of Massey [23]. The 
reaction mixture contained in 1-mL cuvette the appropriate 
amount of the protein, 40 JtM Indigo Carmine or 15 
phenosafranin, 0.4 mM xanthine, 5 |a,M benzylviologen in 
50 mM sodium phosphate buffer with 0.1 mM EDTA (pH 7.0). 
After flushing by argon (99.9999%, (v/v)), a catalytic amount of 
milk xanthine oxidase (0.16 |J,M) was added to the mixture and the 
reduction of FerBHisg-bound FMN and the dye was recorded by 
UV-VIS spectrophotometer (UltraSpec 2000, GE Healthcare). 
For each time, the concentrations of oxidized and reduced FMN 
and dyes at 450 nm and 521 nm (phenosafranin) or 610 nm 
(Indigo Carmine) corrected for the contribution of the appropriate 
dye and FMN, were quantified. 

Rapid Kinetics Measurements 

Stopped-flow spectrophotometr)' was performed on a SFM- 
3000 (Bio-Logic) stopped-flow system equipped with a 0.8- x 0.8- 
mm cuvette (FC-08). The temperature was kept constant at 10 °C 
with a thermostating water bath. Six measurements were averaged 
for each sample. In a typical run 75 |iL of a 30 ^iM protein 
solution was mixed with 75 |J,1 of 10 mM NADH. A flow speed of 
7 ml/ s resulted in a dead time of 0.4 ms. Results from the 
stopped-flow experiments were analyzed in the Bio-Kine software 
(BioLogic). 

Stereospecifity Determination 

[4A-^H]NADH was prepared by mixing of 5.6 mM NAD , 
48 mM deuterated ethanol-de (99,5% D, Sigma), 25 units of 
horse-liver alcohol dehydrogenase, 17 units of yeast aldehyde 
dehydrogenase in 10 mL of 6 mM Bis-Tris/Propane Buffer 
(pH 9.0) at 25 °C. The pH was maintained at constant level of 
the buffer with additions of KOH. The end of the reduction was 
reached during 1-2 hours and was confirmed at 340 nm. The 
enzymes were removed on ice by ultrafiltration through YM3 (cut- 
off limit 3kDa) filters (MUlipore). The temperature of the filtrate 
was decreased to -20 °C and [4A-^H]NADH was lyophilized by 
LyoQuest (Telstar). [4B-^H]NADH was prepared by mixine of 
4.5 mM NAD"^, 5 mM ATP, 6 mM deuterated D-glucose-l-d 
(98% D, Aldrich), 10 mM MgS04, 3.5 units of yeast hexokinase 
and 4 units of glucose-6-phosphate dehydrogenase in 10 mL of 
6 mM Tris-HCl buffer (pH 8.0) at 25 °C. The removing of the 
enzymes and lyophilization were described above. Enzymes and 
nucleotides were purchased from Sigma. The speciphicaUy 
deuterium-labeUed reduced coenzymes were reoxidized by 



FerBHisg. The total reaction volume was 1 mL of 50 mM sodium 
phosphate buffer (pH 7.4) and contained approximately 5 mM 
[4A-^H]NADH or [4B-'^H]NADH, 10 mM 1,4-BQ, (1,4-benzo- 
quinone) and 5 |iM enzyme. When there was no further decrease 
at 340 nm, the mixture was ultrafiltrated and lyophiUzed as was 
described above. The lyophilized samples of NAD"^ were dissolved 
in deuterium oxide (99.8%, Sigma) and aU NMR spectra were 
measured at 25 °C on a Bruker Avance 500 MHz spectrometer 
(Bruker). The control spectra of 10 mM l,4-BQ_and 1,4-BQH2 
were also recorded and no interference was observed in the area of 
the H-4C resonance signal. Chemical shifts are expressed as 
p.p.m. relative to trimethylsilan. 

In-silico Docl<ing 

AutoDock Vina vl.1.2. was used as a docking tool to generate 
putati\'e NAD^/FerBHisg complex. Input files wf;re prepared 
using the Autodock Tools where Gasteiger charges were assigned to a 
Ugand and receptor and subsequentiy non-polar hydrogens were 
merged. The FerBHisg receptor was set as rigid during docking 
and the grid box with dimensions of 15 X 15 X 15A was placed 
into the NAD"*" binding pocket. Exhaustiveness of the search was 
set to the value of 20 while all other AutoDock Vina parameters were 
kept as default. For FMN/FerBHisg complex generation, the 
FMN coordinates were removed from the FerBHisg receptor file. 
The grid box was with dimensions of 40 X 40 X 40 A was placed 
into the FMN binding pocket, the receptor was set as rigid and aU 
other AutoDock Vina parameters were kept as default. 

Results 

Crystal Structure of FerBHise 

The crystal structure of FerBHisg in solution was determined by 
the multiple wavelength anomalous dispersion method using 1 .40 
A resolution diffraction data (Table 1) from crystals of a 
selenomethionine derivative (PDB ID: 3U7R). The crystals belong 
to the orthorhombic space group P2i2i2 with two protein 
molecules in an asymmetric unit. The final model contains 
residues 2-182 (monomer A and B), one polyethylene glycol 
fragment, 451 water molecules and two flavin adenine mononu- 
cleotides. Most residues in the model are well defined in the 2F„ - 
Fc electron density map; missing density at the C terminus 
corresponds to the hexahistidine end (residues 184-190). The two 
monomers are identical with an r.m.s.d. (root mean square 
deviation) between equivalent C-a atoms of 0.08 A. Each 
monomer has an os/pi twisted open-sheet structure containing 
alpha helices on both sides of the central /(-sheet in an 
arrangement typical for a flavodoxin-like fold. Five parallel strands 
in the order (32, (31, P3, P4, [35, are sandwiched between the 
helices al and a5 on one side and the helices a2, a3 and a4 on the 
other side. There are two 3io helices (residues 40-42 and 48-53) 
between P2 and a2 and a short 3io helix (residues 149-151) 
between (35 and a5 (Figure 1). 

Quaternary Structure of FerBHise 

The two molecules of the asymmetric unit form a head-to-tail 
homodimer with two symmetrical active sites facing the opposite 
sides and located at the interface between the- subunits. Dimer- 
ization occurs mainly via anti-parallel packing of helices a3 and 0t4 
into a four-helix bundle. A total of 23 residues from each 
monomer is involved in dimer formation (contact criteria Van der 
Waals overlap > —0.4 A). The interface area of this particular 
assembly covers 11 15.9 A^ according to a surface area analysis by 
PISA [24], which strongly favors it against two other possible 
dimeric assemblies generated in the crystal lattice (interface areas 
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Figure 1. Overall structure of FerBHiSe from P. denitrificans. The two protein subunits are shown in yellow and green ribbon representation 
and two bound FMN cofactors as stick models. The secondary structural elements (five a-helices, three 3,o helices, five (J-sheets and five loops) are 
designated cx, r|, |3 and L and numbered as they appear along the polypeptide chain. 
doi:1 0.1 371 /journal.pone.0096262.g001 



739.2 or 397.1 A^). Moreover, the adjacent dimers might 
possibly interact through flat parts of their surfaces and pack 
perpendicular to one another to form a tetramer (buried area 
12586.5 A^). 

In order to determine the actual oligomeric state of FerBHise in 
solution, size-exclusion chromatography was first carried out on a 
calibrated Superose 12 column. The protein eluted at a volume of 
9.5 mL (Figure SI), corresponding to 1.88 subunits with molecular 
mass of 2 1 288 Da organized in a single molecule. The prevalence 
of a dimeric form of FerB at moderate dilutions agrees with the 
results of earlier studies using chemical crosslinking or light 
scattering [13]. The behavior in more concentrated solutions was 
assessed by SAXS. The theoretical SAXS profiles of the tetramer 
and dimer were calculated from the crystal structure by the 
CRYSOL program and compared with the experimental scatter- 
ing curves obtained at protein concentration of 5, 2.5 and 
1.25 mg/mL (Figure S2A). For concentrations 5 and 2.5 mg/ml, 
the experimental data were better fitted by the tetramer (xr, = 1 .69 
and Xa..-) = 1-27) than the dimer (xs = 8.53 and X2.5 = 5-02) model, 
while for the lowest, both fits gave nearly the same values 
Xi.2.'j= 1-22 for the tetramer and Xi.2.'i= 1-50 for the dimer. This 
diflerence indicates the occurrence of a dimer-tetramer equilibri- 
um. Analysis of the same original data for all concentrations with 
the program OLIGOMER estimated the partial volume fraction 
of dimer:tetramer to be 9:91% iy^ = 1.18) at protein concentration 
5 mg/mL, 13:87% (x^= 1.18) at protein concentration 2.5 mg/ 
mL and 34:66% (x^ = 1.60) at protein concentration 1.25 mg/mL 
(Figure S2B). 

Similarity to Other Proteins 

A search for structurally similar proteins using the structural 
database comparison program DALI returned about 9 hits with Z- 
scores above 15. The top 11 matches belonging to the protein 
family PF03358 and having at least 14% sequence identity with 
FerB are listed in Table S3. Amino acid sequences of these 
matches were aligned with ClustalW and PHYLIP was used to 
draw the phylogenetic tree (Figure 2). FerB clustered with two 
recently structurally characterized proteins, the Escherichia coli 
chromate reductase ChrR [25] and the Gluconacetobacter hansenii 
chromate reductase [26], and also with the Pseudomonas aeruginosa 
FMN reductase [27]. As it is shown in Figure 2 A, all the 12 
representatives included in the analysis share the similar conserved 



sequences thought to be involved in binding of flavin and 
NAD(P)H [27]. 

Flavin Cofactor and its Binding Site 

Contrary to the initial experimental hints [5], the present X-ray 
crystal data show unequivocally that FMN and not FAD is a true 
cofactor of FerB (Figure S3). The isoaUoxazine ring of FMN in 
enzyme crystal is bent 12.8° along the N5-N10 axis. It is uncertain 
at present whether the initially planar oxidized isoaUoxazine 
adopted this conformation as a result of binding to protein or 
whether it was only later converted to a bent semiquinone by 
photoelectrons generated during data collection [28]. The FMN 
binding site is mainly composed of residues from the loops L 1 and 
L3, supplemented by a few residues belonging to the adjacent 
subunit (Tyr46, Asp48, Arg95). The proposed interactions are 
summarized schematically in Figure 3. LI carries a variant of the 
(T/S)ZILQSX(T/S) motif for phosphate binding [27], namely 
SLRKDSLN, of which Serll, Argl3, Serl6, Leul7, and Asnl8 
can form hydrogen bonds with the nonbridging oxygen atoms of 
the FMN phosphate. The amino acid residues 77 to 80 of L3 are 
expected to participate as foUows. Glu77 and Arg80 are bonded 
with one another and maintain two molecules water in positions 
which allow to protonate the 02 atom of isoaUoxazine. 
Alternatively, the proton for 02 can come from the hydroxyl 
group of Serl 13. Tyr78 probably binds FMN through a stacking 
interaction with the dimethylbenzene end of the isoaUoxazine ring; 
the phenolic hydroxyl group of Tyr78 could have a hydrogen 
bond with a phosphate oxygen atom via a molecule of water. The 
peptide backbone nitrogen atom of Asn79 and Arg80 can interact 
with N5 and 04 of isoaUoxazine. As a result, FMN is anchored 
into the pocket in such a way that the pyrimidine end on the 
isoaUoxazine ring remains accessible to the solvent and thus 
available for interacting with electron donors and acceptors. 

Preparation and Characterization of Apoenzyme 

The immobilized metal-affinity chromatography method of 
apoflavoprotein preparation, originally developed for the flavin- 
containing PAS domain of the NifL protein from Azotohacter 
vinelandii [29], proved to be convenient also for FerBHisg, although 
the concentrations of urea and KBr had to be raised to 4 M for 
both in order to remove FMN most effectively. The procedure 
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Table 1. Data collection and refinement statistics. 







3W-MAD data, Se-Met 






Hi-res 


Data set 


Se-peak 


Inflection point 


High-energy remote 


Se-Met 


Space group 


P2,2i2 


P2,2,2 


P2i2i2 


P2,2,2 


Unit-cell Parameters (A) 


0 = 61.22 


0 = 61.25 


0 = 61.25 


0 = 61.03 




f) = 89.19 


b = 89.22 


b = 89.24 


b = 89.14 




c=71.46 


C=71.49 


c= 71.53 


c = 71.28 


Wavelength (A) 


0.97714 


0.97753 


0.97522 


1 .000036 


Resolution range (A) 


30-1.75 (1.78-1.75) 


30-1.75 (1.78-1.75) 


30-1.75 (1.78-1.75) 


19.11-1.40 (1.44-1.40) 


Rsym (%) 


5.7 (33.7) 


6.0 (41.2) 


5.2 (41.7) 


5.5 (28.0) 


//ct/ 


14.8 (2.2) 


14.0 (1.8) 


10.7 (1.4) 


19.5 (3.3) 


Completeness {%) 


86.1 (68.0) 


85.8 (66.0) 


80.5 (57.7) 


97.2 (86.2) 


Redundancy 


3.8 (1.9) 


3.8 (1.9) 


2.2 (1.3) 


5.7 (2.2) 


Wilson B factor (A^) 


21.8 


22.9 


23.7 


15.0 


Refinement 


Resolution (A) 








19.11-1.40 (1.44-1.40) 


No. Reflections 








72133 (4622) 


Rwork/fifree (%) 








16.0 (21.4)/17.4 (22.4) 


No. Atoms 


Protein 








2848 


Ligand/ion 








90 


Water 








451 


B factor (A^) 


Protein 








9.368 


Ligand/ion 








11.465 


Water 








20.158 


RMS deviations 


Bond lengths (A) 








0.013 


Bond angles {^) 








1.765 


ESU (A) 








0.033 



doi:1 0.1 371 /journal.pone.0096262.t001 



yielded apoenzyme with less than 0. 1 % residual enzymatic activity 
and no appreciable absorption in the visible region. At least 99% 
of the original activity could be recovered upon incubation with a 
surplus of the native flavin. Elution volume in the gel permeation 
chromatography (Figure SI) of apoFerBHisg were similar if not 
identical to that of FerBHisg, demonstrating that the apoenzyme 
lacking FMN still retains the ability to form a homodimer. The 
deflavinated product showed decreased storage stability and 
precipitated slowly after being transferred to low ionic strength 
buffers. When subjected to differential scanning calorimetry, it 
gave an irreversible endothermic peak at 56.5 °C, which contrasts 
with a significandy higher temperature of 81.2 °C for denaturation 
of the FMN-containing holoenzyme (Figure S4). Thus, besides its 
catalytic role, the flavin cofactor also displays a structural function 
by stabilizing the dimeric FerB protein. 

Binding Affinity of Apoenzyme for Flavin 

Consistently with the behavior of many other flavoproteins, 
binding of free FMN to apoFerBHise and its mutant forms was 
accompanied by fluorescence quenching, which provided a means 
to determine the ligand dissociation constants (KJ). The Kj values 
for complexes of apoFerBHise with FMN and with riboflavin, 
which lacks the 5 '-phosphate group of FMN, were found to be 



27±2 nM and 6.58±0.05 nM respectively (25 mM Tris-HCl, 
pH 7.4, 25 °C). A 244-fold increase in Kj for riboflavin indicated 
that interaction of apoenzyme with the FMN phosphate signifi- 
cantly stabilizes the complex. This conclusion could be further 
substantiated through mutating amino acid residues, located near 
the phosphate terminus of the bound FMN, to alanine. The 
changes in the Gibbs free energy of binding (AAG4) upon mutation 
were significant for Serll, Argl3, Serl6 and Asnl8 (Figure 4), 
giving quantitative support to the idea that these amino acid 
residues constitute the binding site for phosphate. Similar 
mutational analysis demonstrated the importance of the tyrosine 
residues Tyr46 and Tyr78 and possibly also Serll3 for capturing 
the isoalloxazine ring of FMN. 

Effect of Mutations on Kinetic Parameters for Oxidation 
of NADH 

FerB is known to catalyze oxidation of NAD(P)H by a variety of 
electron acceptors. To identify residues involved in NADH 
recognition and catalytic activity, we determined kinetic param- 
eters of the mutant holoenzymes in the presence of a fixed 
concentration of UQ;0 (0. 1 mM) and variable concentrations of 
NADH. Mutation to alanine of Serll, Argl3, Serl6, Asnl8, 
Tyr46, Tyr78 and Asn79 did not impair the reaction at all. 
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Figure 2. Tertiary structure-based sequence alignment (A) and phylogenetic tree (B) of selected structural homology proteins of 
FerB from a NADPH dependent FWIN reductase family (PF03358). Strictly conserved residues are shown in black, dark gray shades indicate a 
high degree of conservation, while light grays indicate various degrees of low conservation. Generally, shaded residues indicate that a conservation of 
similar residues persists across at least 60% of the alignment. The alignment was generated using ClustalW with a BLOSUIVI62 matrix and default 
parameters and processed with BoxShade 3.31. Organisms and ORF designations from the corresponding genomic sequence: GH-ChrR_G. hanseii 
(D5QFC5, gi: 388603906), ChrR_£ coli (P0AGE6, gi: 84028020), T1501_P. aeruginosa (Q9I4D4, gi: 81541544), FerB_P. denitrificans (A1B9E3, gi: 
69933457), ArsH_S. flexneri (Q7UC03, gi: 75386123), ArsH_S. meliloti (Q92R45, gi: 81634873), WrbA_£ coli (P0A8G6, gi: 67475535), WrbA_P. aeruginosa 
(Q9I509, gi: 81622522), YhdA_S. subtiiis (007529, gi: 81341002), EmoB_Edb (EDTA-degrading bacterium) BNC1 (Q9F9T2, gi: 75466236), Lot6p_S. 
cerevisiae (Q07923, gi: 74583672), and and WrbA_S. pneumoniae (Q97NR6, gi: 81531538). The phylogenetic tree with the greatest bootstrap values 
was constructed by a complete phylogenetic analysis package PHYLIP. 
doi:10.1371/journal.pone.0096262.g002 



whereas mutations at residues Glu77, ArgSO, Arg95, SerllS and 
Glyll5 were effective in decreasing and for most cases in 
increasing K„ (Table 2). The Glu77-Arg80 pair located near the 



pyrimidine side of the isoaUoxazine ring of FMN proved to be 
fundamentally important for enzyme functioning, since elimina- 
tion of either member by charge-reversal or charge-neutralizing 



O 




Figure 3. Predicted interactions between FWIN and the protein. Hydrogen bonds are shown by dashed lines with the bond length (A) printed 
in the middle, while the spoked arcs represent protein residues making nonbonded contacts with the ligand. The image was generated using the 
program LIGPLOT, version 4.5.3. 
doi:1 0.1 371/journal.pone.0096262.g003 
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S11A R13A S16A N18A Y46A E77A Y78A N79A R95A Y97A S113A 

Figure 4. Changes in the standard Gibbs energy for FiVIN binding to alanine substitution mutants of FerBHise. The depicted values of 
AAGf, were calculated from the dissociation constants for the mutant and wild-type protein as determined by fluorimetric titration (see Materials and 
IVIethods). 

doi:1 0.1 371 /journal.pone.0096262.g004 



mutation (Glu77 to Lys, Leu or Met or ArgSO to Glu, Leu or Met) 
reduced 4ai by 2 to 4 orders of magnitude. The observed clianges 
in the apparent values imply that NADH binding was mostly 
affected when mutations occurred at positions Arg80 and Glyl 15. 

Flavin Cofactor Reducibility 

Upon anaerobic exposure to xanthine /xanthine oxidase system 
at pH 7.0 the yellow flavin chromophore of FerBHisg was 
gradually reduced to its hydroquinone form without any 



noticeable formation of a blue or red semiquinone radical 
(Figure 5A). The inclusion of indigo carmine (i?°=— 125mV, 
[30]) enabled us to determine the standard redox potential 
difference between the flavin cofactor and the reference dye (A£"). 
If a true redox equilibrium is attained, a plot of the log([ox]/[red]) 
ratio of the enzyme versus the corresponding log([ox]/[red]) ratio 
of the dye should produce a straight line having a slope of npl 
and an intercept on the vertical axis of -nplS.E°/b'i, where n 
denotes the number of electrons needed for complete reduction of 



Table 2. Kinetic parameters and standard redox potentials of FerBHiSe for the wild-type and mutant variants. 





Km ((iM) 




kcJKm (mM-^ s ^ 


E„ (mV) 


kobs (s ') 


FerBHiSg 


22±2 


258±6 


11727±1100 


-138±2 


147.01 ±0.35 


FerBHiS6-E77A 


180±30 


17.9±3.0 


99 ±23 


-225±10 




FerBHiS6-E77K 


30±8 


0.75 ± 0.05 


25±7 


-239±7 


0.228±0.001 


FerBHiS6-E77L 


250±30 


0.072±0.006 


0.29±0.04 


-229±8 




FerBHiS6-E77M 


420 ±160 


2.5 ±1.0 


6.0±3.3 


-223±4 




FerBHiS6-Y78A 


40±10 


196±13 


4900±1267 






FerBHiS6-R80E 


30±3 


0.275 ± 0.006 


9.2±0.9 


-230±7 


0.608 ±0.003 


FerBHiS6-R80K 


350±80 


9.5 ±1.4 


27±7 


-148 ± 2 




FerBHiSs-RSOiVl 


780 ±280 


0.302±0.084 


0.39±0.18 


-148±6 




FerBHiSe-RSOL 


580±150 


0.142 ± 0.030 


0.24±0.08 






FerBHi56-R95A 


20±2 


146±4 


7292 ±756 


-217±11 




FerBHiS6-R95E 


70 ± 9 


21.1±0.9 


302±41 


-245±1 


45.15±0.09 


FerBHi56-S113A 


60 ±20 


76.7±7.2 


1278±443 


-239±1 


27.88 ±0.03 


FerBHiSs-GllSF 


410±90 


2.2±0.3 


5.4±1.4 




2.607±0.007 


FerBHi56-G115l 


630 ±70 


1.7±0.1 


2.7±0.3 


-149±3 




FerBHiS6-G118F 


85 ± 6 


77±9 


906±124 






FerBHiSg-GllSI 


38±3 


12±1 


316±36 







The apparent Michaelis constant for NADH and the catalytic constant {k^sd were obtained by the initial velocity kinetic analysis at a constant UQ-0 concentration of 
0.1 mM. The standard two-electron redox potential (f ) was determined spectrophotometrically by equilibrating the protein with a redox dye. The first-order rate 
constant for flavin cofactor reduction [kots) was measured upon rapid mixing the enzyme with 5 mM NADH at 10 C under anoxic conditions. Experimental details on 
these measurements are given under Experimental Procedures section. 
doi:1 0.1 371/journal.pone.0096262.t002 
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Figure 5. Spectral changes during the anaerobic reduction of 
FerBHiSg by the xanthine/xanthine oxidase system. The spectra 
shown (from top to bottom) were recorded at 0, 5, 10, 13, 17, 22, 26, 30 
and 36 min for 48 ^M FerBHIse (A) and 0, 5, 10, 15, 20, 24, 27, 30, 33, 37 
and 40 min for 65 |ilVl FerBHIsg with 40 [iM Indigo Carmine. In (B), the 
inset shows the log(ox/red) for the dye vs. that for the enzyme bound 
FMN revealing a slope of 1.08, which is close to the theoretical value of 
1. The standard potential of FMN was calculated from the y-axis 
intercept. 

doi:1 0.1 371/journal.pone.0096262.g005 

flavin (F) or dye (D). Plotting the data in this way (Figure 5B, inset) 
gave a slope of 1, which is close to unity as expected for 
np= nj) = 2. From the intercept value of —0.45 the of bound 
FMN could be calculated as -125 + ((59/2) x -0.45)=- 
1 38 mV. The experiment of Figure 5 was repeated with mutant 
enzymes and, when appropriate, with the alternative reference 
dyes phenosafranin {E°= —252 mV, [30]) or anthraquinone-2- 
sulfonate (ii°= — 225 mV, [31]). The results are listed in the 
penultimate column of Table 2. Testing our system with 50 [iM 
free FMN, we determined an ^"of — 212 mV which is close to the 
published value of —207 mV [32]. The positive shift in E° for 
FerBHisg means that the wild-type apoenzyme binds the oxidized 
flavin more weakly than the reduced form, which facilitates the 
cofactor reduction. This preference is lost in all of the tested 
mutants at Glu77, in the SI 13A mutant and also in the R80E and 
R95E mutants where the positively charged guanidinium group of 
arginine is replaced by the negative carboxyl group of glutamate. 

The reduction of the FMN group by NADH, a plausible 
physiological donor, was evaluated spectrophotometrically in a 



stopped-flow apparatus under anaerobic conditions and at 10 °C. 
Changes in flavin absorbance at 450 nm upon rapid mixing of 
30 pM FerBHisg with 1-5000 1.1M NADH (final concentrations) 
followed single exponential decay kinetics. The fitted value of k„/„ 
responded hyperbolically to the NADH concentration with the 
half maximal effect at 6.4± 1 .0 |a,M and tended toward a limit of 
147.0±0.4 s"' that can be considered as an apparent electron- 
transfer rate constant for the FerB-NADH complex. The 
determined at saturating concentration of NADH (5 mM) did not 
depend on pH in the interval 3-9. By substituting 99% deuterium 
oxid(; for normal water in buffer of pH 6.8, a kinetic solvent 
isotope effect A:„4,(H20)/A:„4,(D20) of 1.9 was noted. The A:,,^, at 
5 mM NADH for the mutant proteins were found to be 
diminished in cases where the enzymatic reaction of NADH with 
quinone was also impaired by the mutation (Table 2). Attempts 
were also made to measure the oxidation of the flavin cofactor, 
pre-reduced with an equivalent of sodium dithionite, by 200 )J,M 
UQjO, but the reaction of all proteins was too rapid and completed 
within the dead time of the stopped-flow instrument (0.4 ms). 

Stereochemistry of Flavin Cofactor Reduction 

Knowledge of the crystal structure of the FerB holoenzyme was 
enough to predict hydrogen transfer to and from the si-ia.ce of the 
bound FMN (Figure 1). However, because our attempts to grow 
crystals of FerB in complex with NADH or NAD"*" had failed, we 
were stiU left uncertain as to the stereospecificity of hydrogen 
removal from C4 of the dihydronicotinamide moiety of NADH. 
Therefore, we examined the stereochemical outcome of the FerB- 
mediated oxidation of NADH that was stereospecificaUy labeled 
with deuterium at either the pro-R or pro-S position. Only in the 
former case the 'H NMR spectrum of the NAD""^ product 
contained a signal at 8.8 p.p.m (marked by arrow in Figure S5A) 
typical for the 4-position of nicotinamide, so we could conclude 
that it is the pro-S hydrogen that is transferred to FMN. 
Quantification of the kinetic isotope effects (KIEs) for the 
enzyme-catalyzed UQ;0 reduction (NADH, 0.15 mM; UQ;0, 
100 nM; FerBHisg, 0.018 nM; pH, 7.4; 30 "C) and for the 
enzyme-bound flavin cofactor reduction (NADH, 1 mM; FerB- 
Hisg, 50 \iM; pH, 7.4; 10 °C) using the C4 monodeuterated 
enantiomers of NADH gave the average KIE ratios of 5.45±0.21 
and6.15±0.34for [4S-%]NADH and 1.07±0.13and 1.05±0.05 
for [4R-^H]NADH. The lower reactivity of the deuterated S- 
enantiomer can be taken as confirmation of the proposed 
stereochemistry. 

FMN and NADH Doclcing 

The laboratory-based studies on interaction of the FerB protein 
with FMN and NADH were complemented with an in silico 
analysis using the AutoDock Vina suite. Initially, in order to assess 
the reliability of the docking method, FMN was picked up from 
the holoenzyme atomic structure (PDB ID: 3U7R) and then 
idealized coordinates of FMN was docked back into the rigidly 
fixed binding site. As illustrated in Figure S6, the docked 
conformation of FMN with the lowest binding energy closely 
resembles its crystal structure, with r.m.s.d. value of 1.1 A for 31 
non-hydrogen FMN atom pairs. The calculated Kj value of 
9.9 nM agrees approximately with that obtained experimentally 
(27 nM, see above). NADH was docked into to the crystal 
structure of holoenzyme where reduced form of FMN was 
replaced by its oxidized form to resemble the proton transfer 
reaction. The complex with the lowest binding energy that 
matches the experimentally determined stereochemistry was 
chosen from approximately 100 docking results. In the putative 
complex (Figure S7) the nicotinamide ring of NADH is stacked 



PLOS ONE I www.plosone.org 



9 



May 2014 I Volume 9 | Issue 5 | e96262 



Structural and Biochemical Study of FerB 




0 -I 1 1 1 1 

0 0.05 0.1 0.15 0.2 

1/[NADH] (^M'^) 



Figure 6. Initial-velocity Icinetics of the NADH:FMN oxidore- 
ductase reaction catalyzed by FerBHise-R95A. Measurements 
were performed at 30 °C with 61 nM enzyme in 25 mM Tris-HCI buffer 
(pH 7.4). Reciprocal initial velocity is plotted against the reciprocals of 
NADH concentration at a series of fixed concentrations of FIVIN equal to 
10 (squares), 20 (triangles), 50 (circles) and 100 (rhombs) |j,M. 
doi:10.1371/journal.pone.0096262.g006 

over the isoalloxazine ring of FMN with a central distance of 4.6 A 
and an angle of their ring planes of 30.0°. The distance of 3.8 A 
between the C4 of NADH and N5 of FMN in this complex allows 
the hydride transfer to occur. The residues ArgSO and Arg95 were 
predicted to contact the nicotinamide moiety. The carboxyamide 
oxygc-n of the NADH (07N) interacts with Arg95 (3.0 A), while 
the oxygen of the ribose (02D) interacts with Arg80 (2.9 A). The 
diphosphate moiety of the NADH could interact via its 02A with 
the amide N of Glyll5 (3.8 A), while the adeninosine moiety 
interacts via 03B with Arg95 (3.2 A) and via N6A and N7A with 
Argl3 (3.4 A resp. 3.3 A). The calculated K,i for the holoenzyme- 
NADH complex (4.4 |iM) is consistent with the experimental K„ 
value for NADH (Table 2). 

Engineering of the Flavin Reductase Activity 

Although FerB belongs to the family of NADPH-dependent 
FMN reductases based on sequence homology, in fact it shows 
very littie activity toward external FMN as compared to that 
toward UQ;0 as a reference, the ratio of the respective kcat values 
being only about 0.006. A possible reason may be a steric 
hindrance preventing the binding of two flavins at once. Inspection 
of the structure suggested that removal of the bulky arginine side 
chain from position 95 might enlarge the active site pocket so that 
the accommodation of the second flavin molecule could be 
facilitated. In line with this expectation, we found that the kcat ratio 
increased to 0.32 as a result of the R95A mutation. The flavin 
reductase activity of the novel mutant enzyme was subjected to a 
bisubstrate kinetic analysis (Figure 6). A double reciprocal plot of 
initial velocities versus concentrations of NADH as the varying 
substrate at four fixed levels of FMN gave a set of parallel lines 
indicating a ping-pong type of kinetic mechanism. By fitting the 
data by non-linear regression to the appropriate kinetic model, the 
catalytic constant and the Michaelis constants for NADH a FMN 
were calculated to be 50±5 s"\ 33±7 ^M and 111±15 ixM. 
Riboflavin and FAD were also substrates exhibiting 4a* values of 
75±10s"' and 38±1 s"' and K„ values of 59±12|xM and 
95 ±7 |j,M. 

Since the side chain of Arg95 protrudes into active site from the 
neighbor subunit across the dimeric interface, more extensive 
changes in quaternary structure had to be considered as an 
alternative explanation for the effect of the mutation. However, 
the persistence of a dimeric state of the mutant proteins R95A and 



R95E, as verified by gel permeation chromatography (Figure SI) 
and SAXS (data not shown), renders this possibility unlikely. 

Discussion 

Proteins of the flavodoxin type display a large range of FMN 
binding affinity. The apoproteins of typical flavodoxins bind FMN 
with a subnanomolar or even picomolar Kj [33], whereas the 
WrbA protein of E. coli has a Kj as high as 2 JiM [34] . FerB 
{Kj=22viM) thus represents an in-between case. The much 
weaker binding of FMN to WrbA compared with that seen for 
flavodoxins was attributed to alterations of critical residues located 
near the plane of the isoalloxazine ring of FMN and to the reduced 
number of hydrogen bonds with the phosphate group [35]. For the 
apoflavodoxin from Anabaena PCC 7119, pivotal binding studies by 
Lostao et al. [22] exploiting not only FMN but also its shortened 
analogues riboflavin and lumiflavin led to a dissection of the 
overall Gibbs binding energy of FMN (—55 kj mol ) into three 
components related to the phosphate (—29 kJ mol~'), isoalloxa- 
zine (—21 kJ mol '), and ribityl (—5 kJ mol ') parts of the 
molecule. An analogous calculation with the experimental 
values for FMN (22 nM) and riboflavin (6.6 \lM) shows that only 
— 13.6 kJ mol ' of the total —43.2 kJ mol ' apparently comes 
from the phosphate, i.e., that the relative importance of 
isoalloxazine may be higher for FerB than for the Anabaena 
flavodoxin. In most flavodoxins isoaUoxazine is sandwiched 
between the aromatic residues tyrosine and tryptophan at the ef- 
face and re-face, respectively [36] . In FerB only the re-face residue 
remains aromatic (Tyr78) and, as indicated by mutant analysis, 
adds to the binding, the contribution of the «-face arginine (Arg95) 
being negligible. Yet another tyrosine residue (Tyr46), which 
extends from the second subunit and has its phenolic ring 
positioned nearly perpendicularly, appears to interact with the 
flavin at its dimethylbenzene end by providing a hydrophobic 
environment. The binding of the flavin ring is further strengthened 
by hydrogen bonds with Serll3 and two backbone NH groups 
(not quantifiable by alanine mutagenesis). Alanine substitutions 
were also helpful in confirming the presumed role of Serl 1, Argl3, 
Serl6 and Asnl8 in binding the phosphate moiety. The sum of 
AAGj values of the individual mutations equals 24.7 kJ mol 
which greatiy exceeds the value of 13.6 kJ mol ' expected from 
the FMN and riboflavin binding data. We interpret this 
discrepancy to mean that the AAGj values cannot be simply 
added because there are some cooperative interactions among the 
spatially close residues. 

Detailed structural studies of NADH binding to FerB-like 
(■nzymes are generally hampered due to the lack of co-crystals 
suitable for X-ray diffraction. At present, the only available data 
are those of the NADH:FMN oxidoreductase (EmoB) from 
Mesorhizobium sp. BNCl [37] that has 21% sequence identity with 
FerB. In the crystal structure of the EmoB-FMN, NADH complex 
(PDB ID: 2VZJ), the ,vf-face (B-side) of the nicotinamide ring of the 
pyridine nucleotide packs against the «-face of the isoaUoxazine 
ring of FMN, as is likely to be the case for FerB. The amine group 
of Lys81 and the backbone carbonyl group of Glyl 12 in EmoB 
were considered to bind to the carboxamide and phosphate groups 
of NADH, respe[:tively [37]. According to the results of amino acid 
sequence comparison (Figure 2 A), the amino acid at position 81 is 
not conserved and EmoB is the only member that has a lysine at 
this site. The equivalent residue in FerB, Asn79, was mutated to 
Ala without affecting enzyme activity. On the contrary, Glyl 12 of 
EmoB is stricdy conser\'ed among all the sequences. The 
corresponding Glyl 15 in FerB is located on a turn between 
strand (34 and helix a4 within a sequence motif SXGXXG that is 
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reminiscent of the common GXGXXG fingerprint region of the 
classical Rossmann fold. These glycine-rich regions are typically 
involved in positioning dinucleotides in a correct conformation by 
allowing for close contact of the polypeptide chain with its central 
diphosphate portion [38]. Our finding that introduction of a large 
phenylalanine or isoleucine substitution at amino acid 1 1 5 results 
in a one-order increase in for NADH and two-order decrease 
in kcat (Table 2) is consistent with (but does not in itself prove) the 
idea of Gly 1 1 5 being a part of the NADH binding site. In our 
model of the FerB-NADH docked complex, Gly 115 lies indeed 
within hydrogen-binding distance of a 02-phosphate oxygen of 
NADH and clear steric hindrances occurs following superimpo- 
sition of Phe or De on it. The same Phe or He substitutions at the 
more distant Gly- 11 8 were much less effective. ArgSO may also 
participate in NADH binding, since its substitutions strongly 
affected the K„, value. In this case the docking model predicts 
interaction with the hydroxy! group of NADH ribose. 

New structural data link FerB very closely to the recently 
characterized chromate reductases (ChrR) of Escherichia coli (PDB 
ID: 3SVL) and Gluconacetobacter hansmii (PDB ID: 3S2Y). One thus 
would expect that all three enzymes might operate by the same 
overall mechanism, but this need not necessarily be true. 
Previously it was claimed that the reduction of chromate by 
NADH catalyzed by G. hansenii ChrR proceeds via an ordered, 
sequential reaction with chromate binding first followed by 
NADH and that binding of NADH alone produces a death-end 
complex [26]. Some facts about FerB, however, contradict this 
picture. When catalyzing a reaction of l)c't\vecn NADH and 
quinone, FerB obeys a ping-pong type reaction mechanism [5], 
which is consistent with a temporary parking of reducing 
equivalents from NADH at the enzyme-bound flavin cofactor. 
Prestationary assays reported here show that the flavin group of 
FerB is still reducible by NADH despite the absence of any 
electron acceptor like quinone or chromate. The estimated 
limiting k„i,s (~140 s~' at 10 °C) is similar to the steady-state 
value of A;^^, for the NADH-driven reduction of UQ;0. Taking also 
into account a comparable deuterium kinetic isotope effect upon 
both processes, we conclude that (i) an electron-acceptor- 
independent reduction of flavin is operative, (ii) it is kinetically 
competent to account for the observed enzyme activity, and (iii) it 
may be rate-limiting for the whole catalytic cycle, at least with 
UQ;0 as electron acceptor. 

Following hydride transfer to the N5 position of the flavin, 
stabilization of the negative charge evolving at the Nl-C2 = 02 
locus of the isoalloxazine moiety is required. In flavoenzymes this 
is usually accomplished through the interaction with a positively 
charged entity, either fully charged, such as a Lys or Arg side 
chain, or partially charged, such as the N terminus of a heUx or a 
cluster of peptide nitrogen [39] . In the canonical crystal structure 
of the rat liver quinone reductase (PDB ID: IQRD) 02 of the 
FAD cofactor is hydrogen-bound to the bac:kbone amide groups of 
Glyl49 and GlyL5() and also to the side chain of Tyrl55. It was 
postulated that Tyrl55 and the neighboring His 161 form a 
charge-relay system needed for the reaction to take place without 
unfavorable charge separation [40] . FerB and the ChrR proteins 
of E. coli and G. hansenii are distinguished by a different 
microenvironment of the Nl-C2 = 02 locus, comprising a 
hydroxyl group of serine and a water molecule hold in position 
by a glutamate-arginine pair. The residues Sc-rl 13 and Glu77 of 
FerB are well conserved among the homologous sequences 
(Figure 2A).The importance of Serll3 gains support from 
experiments of the alanine mutant that exhibits lower specific 
activity as well as depressed FMN binding. The involvement of 
water in catalysis by FerB is consistent with the observed solvent 



isotope effect of 1 .9 for FMN reduction, suggesting that an oxygen- 
hydrogen bond is broken in a rate limiting step. Water has already 
been proposed to serve as a proton donor to Nl or 02 atoms at 
flavin reduction in several other flavoenzymes, such as NADPH- 
cytochrome P4,50 oxidoreductase (PDB ID: lAMO; [41]), 2,4- 
dienoyl-CoA reductase (PDB ID: 1PS9; [42]), FMN-dependent 
azoreductase (PDB ID: 2HPV; [43]) and peptidyl-cysteine 
decarboxylase (PDB ID: 1G5Q; [44]). A crucial role of Glu77, 
amply substantiated by the results of mutational experiments, 
probably relates to maintaining the architecture of the active site. 
Unlike FerB, where a pair Glu77-Arg80 exists, the homologous 
glutamate in both of the ChrR enzymes pairs with a more remote 
arginine located 43 amino acid residues away. In E. coli ChrR this 
arginine lies at the surface and participates in dimer-dimer 
interaction [25]. Subtle differences in hydrogen bond networks in 
ChrR and FerB thus may explain why the former enzyme forms a 
stable homotetramer while the latter remains dimeric in dilute 
solutions. 

Of the proteins similar to FerB listed in Table S3, five have been 
demonstrated to catalytically reduce flavin, i.e., to act as an 
NADH:flavin oxidoreductase. Judging from their crystal struc- 
tures, these enzymes have the flavin cofactor largely exposed to the 
bulk solvent, whereas in FerB the dimethylbenzene part of the 
flavin nucleus appears to be shielded by the side chain of the 
Arg95 residue from the neighbor protomer within the homodi- 
meric structure. The mutant protein R95A, prepared from 
FerBHisg to reheve the presumed steric constraint, indeed gained 
the ability to reduce FMN with a kcat of 50 s~\ which is close to 
the value for the B. subtilis YhdA (60-70 s"', [45]) but lower than 
that for EmoB of Mesorhizobium sp. bacterium BNCl (397 s~', 
calculated from the V^^^x of 1 1 1 1 pmol min ' mg ' given in [37], 
using a subunit molecular mass of 21 412 Da). Similarly to EmoB, 
the engineered FerB obeyed a ping-pong kinetic mechanism with 
NADH and FMN as substrates (Figure 6) suggesting that NADH 
and a loosely associated FMN molecule interact in an alternating 
fashion with the same pocket containing another, tightly bound 
FMN. Such a mechanism was previously identified for EmoB 
through preparation and structural analysis of its FMN-NADH 
and FMN-FMN complexes [37]. Since a single mutation 
increasing the accessibility of the bound FMN sufficed for FerB 
to become an active flavin reductase with affinities for FMN, FAD 
and riboflavin (indicated by K„ values) comparable each other, 
isoalloxazine rings contacts may be of importance in recognition of 
the external flavins. According to the H NMR chemical shift 
measurements, FMN self-associates in aqueous solutions with an 
association constant in the range 20-100 M ' [46]. This would 
give an expected dissociation constant of the order of 10 mM, 
which is however still two orders of magnitude al)ove the i'stimat(;d 
Michaelis constants (~0. 1 mM). Pre-existing protein environment 
around the bound flavin cofactor thus probably also contributes to 
the external flavin substrate binding. A close relationship of FerB 
to the flavin reductases, demonstrated now both structurally and 
biochemically, may shed some light on the evolution of flavin 
reduction and various promiscuous catalytic activities of flavopro- 
teins in general. 

Supporting Information 

Table SI The mutagenic oligonucleotide primers for site- 
directed mutagenesis of FerB Hisg. 
(DOC) 

Table S2 Data collection fi-om SAXS. 
(DOC) 
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Table S3 Amino-acid identity and structural similarity across 
the structural homologs of FerB in a FMN reductase protein 
family (PF03358). 

(DOCX) 

Figure SI Superose 12 analytical size exclusion chromatography 
of FerBHise holoenzyme (1), apoenz)'me (2), and FerBHis6-R95E 
mutant (3). The calibration curve used to estimate the native 
molecular weight based on the elution position is indicated. 
(TIF) 

Figure S2 SAXS scattering curves. X-ray scattering of 
FerBHisg was compared with calculated scattering of the dimer 
{dashed line) and the tetramer (solid line) computed by CRYSOL for 
concentrations 1.25 (A), 2.5 (B) and 5.0 (C) mg/mL. The program 
OLIGOMER was used to estimate dimer and tetramer volume 
fractions for concentrations 1.25 (D), 2.5 (E) and 5.0 (F) mg/ml. 
(TIF) 

Figure S3 Electron density map of FMN. Stereo view of the 
2Fg - Fc electron density map contoured at 1.0 a cut-off 
corresponding to FMN at the active site. 
(TIF) 

Figure S4 Thermal stability of apoFerBHisg (A) and 
FerBHise (B) followed by DSC. Protein concentrations were 
adjusted to 4.3 mg/mL for apoFerBHisg and 2.5 mg/mL for 
FerBHisg in buffer C and data were collected from 25 to 95 °C at a 
heating rate of 60 °C per hour. 
(TIF) 

Figure S5 NMR spectra of NAD+. The 'H NMR aromatic 
region of NAD+ formed from [4B-^H]NADH (A) or 
[4A-^H]NADH by FerBHisg in the presence of 1,4-benzoquinone 
showing a resonance signal at 8.8-8.9 p.p.m. (A) or no resonance 
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signal (B). The enzyme removed the proton from B-side (pro-6', si- 
face) of the dihydrogennicotinamide ring. 

(TIF) 

Figure S6 Superimposition of FMN from the crystal structure of 
FerBHisg (white) and from the docking-derived complex (yellow). 

(TIF) 

Figure S7 Structure of the docked FerBHise-NADH 

complex. NADH was modeled into the FerBHisg structure using 
AutoDock Vina vl.1.2. FMN and NADH are shown with yellow and 
white carbons, respectively (A). The nicotinamide ring was stacked 
on top of the isoalloxazine ring at distance 3.8 A that 
corresponding with distances of other NADH-dependent FMN 
reductases. The residues are depicted in green and violet to 
represent belonging to two different subunits. The known crystal 
structure of EmoB complex with NADH (PDB ID: 2VZJ) is 
included for sake of comparison (B). 
(TIF) 
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